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ABSTRACT 


A plasma was produced by the interaction of a 300 MW, 25 nsec, Nd 
laser pulse with an aluminum target. The resulting plasm expanded into 
an ambient background of 2.5 X 1072 torr and was analyzed using floating 
double probes, magnetic probes and a quartz pressure probe. An early 
disturbance was noted before arrival of the main plasma. 

Further experiments separate this early signal into photoelectric 
response and two fast plasma pulses traveling with constant speeds of 
bed, X 10° em/sec and 5.9 X 10! emysec. 

Mapping of the plasma density indicates that the early plasma pulse 
is not symmetric with respect to the target normal but expands along a 
line defining the reflected laser pulse. This same mapping indicates 
that the main plasma expands anisotropically for the first 60 nsec 
resulting in an early tine asymmetry with respect to the target normal. 
At times greater than 120 nsec, the asymretry of the main plasma is no 
longer evident and its density distribution is za with ere 


to the target normal. 
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I. INTRODUCTION 


Karly investigations anwar ne laser induced electron emission by 
| William I. Linlor [Ref. 36] in 1962 were followed by a multitude of 
papers dealing with laser produced plasmas. Lichtman and Ready [Ref. 
34] were among the first to detect an early signal preceding the main 
plasma body. They ascribed the signal to reverse photoelectric effect 
[Ref. 34] from the electrostatic probes being used. This hypothesis 
received immediate criticism from Honig (Ref. 23] who ascribed the 
early signal to a fast thermal effect. Since then, several new 
hypotheses including ionizing potential waves [Refs. 2 and 31] and 
detonation or blast waves [Refs. 9 and 22] have been forwarded to 
explain this early disturbance. However, E review of the literature 
does not reveal a conclusive report on the behavior or content of this 
signal. 

The separate discovery of spontaneously generated magnetic fields 
by Korobkin and Serov [Ref. 29] and by Stamper [Ref. 43] further 
stimulated interest in this area. Computer simulations have indicated 
the presence of megagauss magnetic fields early in the plasma expansion. 
Fields of this magnitude will influence plasma containment schemes 
generated in the controlled fusion effort. 

In response to an obvious need for complete understanding of the 
laser produced plasma and its induced magnetic fields McKee [Ref. 38] 
and Bird [Ref. 8] undertook extensive surveys of the self-induced 
magnetic fields and attempted to relate them to the expanding plasma. 


During Bird's investigation he noted the presence of an early magnetic 
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Siegal. This revelation ted to renewed interest in the early electro- 
static probe signals. There are four possible sources for this early 
signal. They are: 

1. laser induced electron emission. 

2. Fast plasma blowoff. 

3. Electromagnetic disturbances. 


4, Ionization of the background gas by plasm produced photons or 
ionizing potential waves. 


This investigation has been conducted in an effort to determine which 
@teemese Causes 1S responsible for the observed signal and to provide 
a more detailed understanding of the behavior of the phenomenon producing 
the signal. 

The thesis is divided into six more sections and two appendixes. 
Seetion Ir is a description of the general experimental arrangement. 
Section III describes the experiments conducted to solve the problem 
and discusses the results. Section IV attempts to relate two of the 
observations made in Section III to theory. Section V presents a 
Ana with conclusions. Section VI provides recommendations for 
further investigation. Appendixes A and B describe the construction 


and theory behind two of the diagnostic tools used during this research. 


e 





II. EXPERIMENTAL ARRANGEMENT 


A. LASER 

The laser used was a (ened K-1500, Q-switched neodymium-in-glass 
system. The system consists of an oscillator which generates the 
initial 25 nsec (full width at half power) pulse through Q-switching by 
Pockels Cell. The oscillator pulse is expanded and passes through an 
amplifier wnich increases the energy of the pulse. The arrangement is 


shown in Figure 1. A detailed description is provided by Davis [Ref. 7]. 


CW Alignment Laser ) 
Pockels Cell a 


Polarizer —(/) 


Kel Nd Oscillator —LY 


Beam Expansion Y 
Optics (.5 to .75 inches) 


A aii ode 
N 


N 
K-2 Amplifier— = ergo Diffuser 
Beam Br. eg E 
/ = 
/ £| [Energ 
‘Focusing g 
Probe SA Lens E ere er Tektronix 
Vacuum Chamber £ 5648 
m Storage 
Tarzet (OM Oscilloscope 
E OR eter Delay Tek. 7704 
Line Oscilloscope 


Figure 1. Block Diagram of the General Experimental Layout 
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The output of the laser was controlled by varying the voltage 
applied to the K-1 and K-2 flashlamps. The range of energies 
obtainable is 2.5 to le Joules. The duration of the pulse is 25 nsec 
giving a power range of 100 to 480 MW. The laser focal spot size is 
approximately 0.15 cm in diameter after focusing. This gives a eee 


(5 Joules) power density of 2.8 X 107 watts/cmÉ. > Laja 


TS 








B. LASER MONITORING TECHNIQUES 

Approximately one percent of the laser pulse energy was reflected 
by a beam splitter onto an MgO block. ‘The laser output was monitored 
from this block with a Korad K-Dl photodiode after passing through a 
0.1 percent neutral density filter. This photodiode provides one signal 
proportional to the pulse power and a second, integrated, signal 
proportional to the laser pulse energy. 

As shown in Figure 1, the power signal was used to trigger a 
Tektronix 7704 oscilloscope and was also displayed on that oscilloscope 
through the positive input in a differential amplifier. The energy 
signal was displayed on a Tektronix 564B storage oscilloscope. 

The photodiode energy signal was calibrated by L. L. McKee [Ref. 38] 
using a Westinghouse RN-1 Laser Radiometer. The radiometer was used to 
provide an absolute measure of the incident pulse energy. 

The energy output of the laser was monitored on every shot and the 
power was displayed on the same oscilloscope used to take probe response 
data. For successive shots during the mapping procedure the energy was 
allowed to vary by + 10%. Shots with energy variations exceeding this 
were excluded from the data base. When properly aligned the laser 


produced very consistent pulses. It should be noted that if shots were 
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repeated at intervals of less than one minute, the laser energy tended 
to increase slightly from shot to shot. A one minute cooling period 


between shots resulted in excellent reproducibility. 


C. TEMPORAL REFERENCE FRAME 

The temporal relation between the diagnostic probe response and the 
leading edge of the laser pulse was determined by taking into account 
the various cable delays and laser beam optical path lengths. The 
geometry is depicted in Figure 1. A long (10 meter) delay line was 
used from the diagnostic probe to the oscilloscope in order to separate 
the laser pulse input and early diagnostic responses which would other- 
wise be superimposed by the differential amplifier. The computed delay 
was 33 nsec assuming the signal to travel at 1.98 X 1019 em/sec in the 
RG-174/U coaxial cable. 

To experimentally verify this delay, the 100 volt calibration signal 
from a Tektronix 555 Dual Beam Oscilloscope was fed into a Tektronix 
Tunnel Diode Pulser. The pulser output was fed into a "T." One side 
of the T went to the positive side of a Tektronix 7A-13 differential 
vertical plug-in on a Tektronix 7704 oscilloscope through the photo- 
diode power signal cable. The other side of the T went to the negative 
side of the same amplifier through the delay cable. A reproduction of 


the oscilloscopic trace obtained in this manner is given in Figure 2. 


A TA 


a 
37 nsec delay 





Fleur 2. Ten meter cable time delay measurement. The horizontal scale 
is 20 nsec per division. The measured delay is 3/ nsec. 
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When the measured 37 nsec delay is corrected for the difference in 

optical path length from the beam splitter to the target and the MgO 
diffuser block the delay is found to be 35 nsec. This experimentally 
measured delay agrees very well with that calculated earlier and is 

the delay used throughout this thesis in correlating diagnostic probe 
response with the laser pulse. Zero time is arbitrarily established 
as that time at which the leading edge of the laser pulse strikes the 


target. This is depicted in Figure 3. 





eee Pulse 


Defined Zero Time 
Figure 3. Definition of Zero Time. The horizontal scale is 20 nsec 


per division. i 
D. PLASMA CHAMBER 
The target was located inside a chamber specifically designed by 
McKee [Ref. 38] to facilitate diagnosis of the laser plasma by probes or 
optical means. A horizontal cross section of this chamber is depicted 


in Figure 4. 
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Target Normal 
(Z-Axis) 


Incident 
Laser Beam 


Aluminum 





Chamber Height 
is Six Inches Scale 


Target Rotation Knob q" 


Figure 4. Top view of vacuum chamber. Port #1 is the laser beam entry 
port, #2 is the reflected laser beam observation port, #3 is 
the transmitted laser beam observation port, ports numbered 
4, 5, 6, and 7 are optical/probe observation ports. There is 
also an optical/probe port in the top of the chamber. 
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The laser pulse enters the chamber after passing through a 
converging lens and strikes the target at an angle of 30 degrees with 
the target normal. This thirty-degree incident angle permits probing 
along the target normal which appears to be the axis of greatest 
symmetry in relation to the expanding plasma. The converging lens was 
normally defocused to minimize cratering of the target. 

The target is a flat disk of one-half inch 6061 aluminum alloy 
plate. It is two inches in diameter. The target was rotated every 25 
wens) asMets tomprevent severe cratering. 

The chamber was pumped to a vacuum of 2.5 X 1972 ere (air) by an 
ol diffusion pump. This pressure was measured with an ionization 
gauge attached directly to the chamber. Unless otherwise stated, all 


measurements were made at this chamber pressure. 


| 
| A Angular Reference Line 
| i 
1@=0 
Target Material 


Laser Focal 
Spot 


Figure 5. Cylindrical polar coordinate system used in this thesis. 





The chamber coordinate system is polar with its center at the focal 
Spot of the laser pulse. This is illustrated in Figure 5. Throughout 
the text reference positions refer to cylindrical polar coordinates. 
Distances will be in centimeters unless otherwise specified. As an 
example, (0.5, 7/2, 1.0) refers to the position r = 0.5 cm, @ = 7/2 and 


z = 1.0 cm. 


E. MAGNETIC PROBES 

Probe construction and calibration is covered by McLaughlin [Ref. 
39]. During this work a single magnetic probe consisting of five turns 
of #40 Formvar copper wire with a diameter of one millimeter was used. 
Using the procedures outlined by McLaughlin, an effective area of 
nA = 3.1 X ieee má was measured for the coil. The coil proved to have 
a linear frequency response up to 12 MHz. 

Since the probe was used only to determine temporal relations and 
to verify azimuthal symmetry of the self-induced magnetic fields, no 
integrator was used and the data was recorded simply as the voltage 


induced in the coil. 


F. FLOATING DOUBLE PROBES 

Appendix A contains details of the double probe construction and 
theory. The primary use of these probes was to determine the time at 
which disturbances arrived at the probe and to determine the relative 
density of the plasma at a given time in a given position. Probes were 
biased sufficiently negative to collect saturation ion currents. 

Probes used in this investigation have been numbered one through 
four. Probe number three has larger electrodes and was used 


extensively due to its higher signal to noise ratio. The characteristic 
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for the peak density of the main plasma pulse is reproduced as Figure 6. 
The characteristic was measured in a background of 2.5 X 10°? torr with 
mom nel laser energy of 5 + 0.5 joules. The characteristic is 
Symmetric and appears to partially saturate at + 10 volts. However, 
there is another, smaller, break in the curve at 20 to 25 volts. 
Therefore, the probe was biased at -20 volts to collect a saturated ion 
current. A large secondary emission coefficient of «= 0.5 was used in 
obtaining a quantitative idea of plasma density to compensate for this 
large biasing voltage. 

By carefully isolating the entire system from ground, recorded 
probe response voltages at zero bias were reduced to less than 300 mV. 
For most probes the response was zero at zero bias indicating that the 
probe was in fact floating with the plasma potential. 

To obtain an estimate of the plasma density, equation (29) of 
Appendix A is solved for n; and evaluated with Ze = 1.5, €= 0.5, 


R= 3.60, A= 4.6 X 10-3 em and 1 given a variety of significant 
stream speeds: 


(1) n (cm2) = 1.59 x10% y (volts) 
i Vip (cm/sec) 
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Figure 6. Probe Number 3 Characteristic at the Peak of the Main Plasma 
Pulse. 
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Table (1) gives n, versus probe response voltage for a number of 


i 
significant plasma stream speeds. 


Vip (cm/sec) a (emi) probe Response (volts) 
1.0 X 10! 1.59 x 1023 
5.0 x 10! 3.18 X 1012 
1.0 X 108 9.75 x 101 


Table (1) Plasma Density Versus Probe Response Voltage for Probe 
Number Three 

G. PRESSURE PROBES 

Pressure probes are useful in measuring the momentum per unit area 
transported by the plasma. In this investigation use was made of the 
piezoelectric properties of an X-cut quartz crystal. The momentum 
transfer caused by plasma impact with such a crystal creates a stress 
wave. Propagation of this stress wave through the crystal results in 
Eninduced voltage. Theoretical details describing this induced 
voltage are included in Appendix B. 

The pressure probe used in this investigation was constructed ard 


provided by the Stanford Research Institute. 


H. DATA COLLECTION ERROR ESTIMATES 
l. Probe Positioning Errors 
The diagnostic probes were inserted into a brass collet system 
which was welded onto a brass base plate. The base plate could be 
moved about on an "0" ring seal located on the top of the chamber with- 


out breaking the vacuum seal. The base plate was scribed with 


reference lines. These lines allowed positioning of the plate over a 
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one millimeter grid attached to the top of the chamber. A new origin 
was determined at the start of each data session and marked on the grid. 
The probe was then moved to new positions using the grid as a reference. 
The position of the probe relative to the target was measured internally 
at the beginning of each data session and again at the end of the 
session to insure that the probe had not moved in its holder. While not 
very sophisticated, this method did allow positioning of the probe to 
within + 1/2 mm of the stated position in relation to the target. 
2. Absolute Plasma Density Errors 

The relative measure of eer density from point to point is 
considered good due to the high reproducibility of the laser and plasma 
pulses. It is estimated that these relative values are accurate to with- 
ay + 10 percent. 

Absolute plasma density calculations are subject to several 
errors. 

a. High Probe Bias 

A double probe bias higher than that at which initial 

saturation occurs may have generated errors as large as 100 percent in 
the recorded probe response due to secondary emissions and plasma 
perturbations. This error is consistent for all probe responses and 
can be ignored in the relative density plots. To help reduce the 
possibility of introducing a significant error in the absolute density 
determinations, a large secondary emission coefficient of €= 0.5 was 


used. 
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b. Deviations In the Laser Pulse Energy 
While the laser was very reliable and gave very reproducible 
pulses, the output energy was allowed to vary by as much as + SIURpercent. 
The early signal probe response to the laser's energy content was found 
to vary linearly. The Slope of the line describing this relation is 
Such that the allowed deviation in laser energy could result in a 
deviation of + 7 percent in the amplitude of the early signal response. 

c. Error in Probe Tip Measurement 

The area of the probe tips is considered accurate to within 
ten percent. | 

d. The plasma stream speeds are considered accurate to within 
sens percent, 

In view of these many deviations plus subtle factors [Ref. 48] 
in interpretation of double probe signals the absolute plasma densities 
are not considered accurate to better than a factor OÍ two. 

3. Errors In Temporal Resolution 

The time resolution in these experiments is considered very good. 
At this background density, the plasma expansion speed is nearly 
constant within the allowed range of energies. The temporal resolution 
on individual photographs is enhanced by inclusion of the laser pulse 
which serves as a positive indication of zero time. Temporal 
correlations and measurements are made with a horizontal scale of less 
than 40 nsec per cm. This allows interpretation of the results to with- 
in 4 nsec. For two double probes 1.75 cm apart, 4 nsec represents 46 
percent of the transit time for the fast pulses travelling at 2 X 10° 
cm/sec and represents only Y percent of the time involved in passage of 


the main plasma front at 2 X 10! cm/sec. When probe spacing and laser 
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ERC rev errors are considered, the fast signal speeds should be good to 
within 47 percent and the main signal speeds to within six percent for 
single data points. Im each case the speed has been averaged from more 
than 25 data points so that the respective limits are + 10 percent for 


the early signal speed and + 2 percent for the main plasma speed. 
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Ill. INDIVIDUAL EXPERIMENTS AND THEIR INTERPRETATION 


The thesis problem was attacked through a series of experiments 


designed to reveal the behavior and make-up of the early disturbance. 


A. DISTURBANCE SPEED MEASUREMENTS 

The first set of experiments was made in an effort to accurately 
determine the speed of component parts of the laser produced plasma. 
These speed measurements were accomplished in two different ways. 

l. Pulse Arrival Time at a Geographical Position 

in this experiment, a series of oscilloscopic photographs were 
taken at points along the z-axis. The series starts at 0.5 cm and 
Sons SES tor data points at one millimeter intervals to 3.0 cm, at two 
millimeter intervals from 3.0 to 4.0 centimeters and a final data point 
at 5 cm. This series was repeated three times. 

a. The first series was taken with the oscilloscope set at 100 
nsec and 2 volts per cm. These settings are designed to enhance 
determination of the speed of the main plasma pulse. A typical probe 
response is traced as Figure 7 with the pertinent features labeled. 
Main plasma front passage times were measured for all data points. 
These times are plotted against the probe position in Figure 8. The 
Slope of the line connecting the data points represents the speed of the 
main plasma front at any point. As can be seen the velocity is nearly 


consistent at 1.1 X 101 cem See within the limits of this survey. 
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Figure 7. Double Probe Response at (0, 0, 1.0). The dashed line 


indicates the time at which the main plasma front passed 
the probe position. 
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Figure 8. Time of Passage of the Main Plasma Front Versus Probe 
Position: 
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Assuming a constant speed and extrapolating backward indicates that 
the main pulse left the target very early after laser pulse impact 
(approximately 10 nsec). 

b. The second series was taken with oscilloscope scales of 
20 nsec and 500 mV per cm. This scale enhances direct measurement of 
the early disturbances. A typical oscilloscopic trate is copied as 


Figure 9. 
Main Plasma Front 
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Figure 9. Early Laser Plasma Produced Disturbances Recorded by a 
Double Probe at (0, 0, 1.25). 


miiewexacy front of Lhe early signal is often difficult to fix. 
This results in the large data scatter seen in early disturbance data. 

A plot of early disturbance front passage versus probe position is 
reproduced as Figure 10. The slope of this plot indicates a fairly 
constant speed of 1.1 X 10° em/sec. Also included in Figure 10 is a 
similar plot for the second early disturbance. The slope of this line 
indicates a front speed of 5.9 X 10! cm/sec. 

It should be noted that the speed of these early disturbances is 
reasonably constant even to the point where the signals are lost in 
system noise. 

By extrapolating backward assuming a constant speed we find that 


the various disturbances left the target at the following times; main 
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plasma at t = 0 to 10 nsec, first early disturbance at about 3 nsec and 


the second early disturbance at t = 12 nsec. 
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Figure 10. Plot of the Time of Passage of the Early Disturbance Versus 
Probe Fosition 


2. Two Double Probes 

To obtain more accurate speed information, a two probe holder 
was constructed to allow simultaneous sampling of the same plasma by 
two tandem probes. The holder is illustrated in Figure ll. 

The probes used were very fine to minimize perturbation of the 
streaming plasma. The probes were placed in the holder 1.75 cm apart. 
They were positioned in the chamber so that one probe was one millimeter 
to the right of the target normal and the other was one millimeter to 
the left of the target normal. Each probe was biased at -20 volts and 


> 


the chamber pressure was reduced to 1.8 X 10 “ torr. A series of 32 


À 


pictures were taken on a Tektronix 555 Dual Beam Oscilloscope. The 


ny I 
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Figure 11. Device for the Insertion of Iwo Double Probes 


scope was set at 500 mV and 40 nsec per cm. A typical oscilloscopic 
trace of the probe response is depicted in Figure 12. At ten of the 
data points, the probe cables were reversed. In each instance, the 
same separation between arrival times was noted. The reversed cable 


trace is also included in Figure le. 
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Figure 12. Simultaneously recorded tracings from two double probes 
sampling the same plasma. (a) The first probe is at 
(0.1, 7/2, 1.1) and the second probe is at (0.1, 37/2, 
2.85). (b) Taken at the same data points but with the 
cables reversed at the oscilloscope. 

The average speed of the early signal as computed from this data 
base is 1.0 X 10° cm/sec. The average speed of the main plasma signal 
momneund to be 1.1 X 10! cm/sec. These results agree very well with 


those found by the first method. A summary of these speeds is given 


in Table 2. 
Plasma Feature speed 
First Early Disturbance 1. ek 10° em/sec 
Second Early Disturbance 5.9. X 10! cm/sec 
Main Plasma Front 1.1 X 10/ cm/sec 


Table 2. Summary of Significant Speeds 
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B. EFFECT OF A TRANSVERSE ELECTRIC FIELD 

A variable transverse electric field of up to 13,000 volts/cm was 
applied to a collimated plasma to determine whether the field would 
deflect the plasma beam. This will reveal information on the charge 
to mass ratio of the disturbance if it is particulate in nature. 

Figure 13 depicts the device constructed to provide a collimting 
effect on the laser produced plasma. The devise was designed with 
several purposes in mind: 

1. To collimate the plasma. 


2. To help shield the electrodes from the plasma and thereby 
reduce the possibility of discharge across the conductive plasma. 


3. To reduce the total number of plasma particles subject to the 
field in order to reduce the effects of shielding the plasma. 


(a) (b) (c) 


62) 


target 
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.32cm pinholes 


Laser Pulse | Plasma 


Figure 13. Collimating Device (a) General View, (b) Top View, and 
(c) Front View 
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Figure 14 depicts the probe response at (0, 0, 5.8) with and without 
the plexiglass collimating device in place. As can be seen, the probe 
response has been altered; however the early disturbance and main plasma 
pulses remain distinguishable features of the response and therefore the 


pinhole collimating device is usable. 


100 nsec/soomV ' 





Figure 14. Effect on the plasma produced by the presence of the 
collimating device. (a) Double probe response without the 
device. (b) Double probe response with the device in place. 
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Figure 15. Electrodes and Their Mounting 
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Figure 15 depicts the electrodes used in this experiment and the 
mounting equipment. All surfaces of the entire structure are heavily 
coated with General Electric "RIV" to eliminate corona formation. 

Assuming that the electrodes form a pair of infinite parallel 
plates (ignore fringing), the stored charge on these electrodes is found 
CO De: 

(2) a=EA Where (A) is the plate area 

(3) V= Ed Where (d) is the electrode spacing 

(4) q = VA | 
With V/d = 13,000 volts/cm and (A) = one cm“ the stored charge is 
ley X 1077 coulombs or TEEN. 10? electrons. 

Plasma densities greater than this will be able to completely shield 
themselves. Typical electron densities in the main plasma pulse are 


approximately 101" cm while probe response to the early disturbance 


ol ion charges/cm. 


indicates a density of approximately 2.4 X 1 
Assuming a nearly neutral plasma in each case it is obvious that both 

pulses are capable of completely shielding themselves. The collimating 
device has pinholes 3.2 millimeters in diameter. This reduces the early 


pulse charge density to 2.5 X 10 om >. 


The early plasma is still 
capable of shielding itself. However, in doing so, nearly all of the 
particles would be lost from the beam resulting in a significantly 
reduced signal. This argument assumes that the plasma particles have 
time to respond to the applied field. 


Because of the shielding problem a large capacitor was used in the 


electrode charging circuit. This circuit is shown in Figure 16. 
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Figure 16. Electrode Charging Circuit 


The RC time for this circuit is approximately 2 nsec between the 
storage capacitor and the electrodes. The capacitor stores 4.5 X 1073 
coulombs of charge at 18,000 volts. Therefore the electrodes will 
receive a significant amount of additional charge from the power supply 
semi che 27 nsec in which the early disturbance is traversing the field. 
dine charge available in the storage capacitor should be sufficient to 
overcome any shielding in the early signal. 

Deflection of singly ionized aluminum atoms assuming no shielding 


of the field by the early disturbance will be: 


(5) D atº/2 Where (D) is the lateral displacement 


(6) a 


F/m = eE/m = eV/md 
If V/d = 13,000 volts/cm then the acceleration is 4.6 X 104 cm/sec? 
for aluminum ions and 8.5 X Lot see for electrons. 
The early disturbance travels with a speed of approximately 
1219 X 10° cm/sec. Therefore, it will be within the 3 cm field for 
t = 27 nsec. Substituting these values into equation (5) we find that 
the ions incur a transverse displacement of 0.9 mm during transit 
while the electron lateral displacement is 170 cm (ignoring space charge 
effects). 
These values indicate that if the early disturbance was composed 


entirely of electrons, it should be completely deflected and no signal 


> 





would be received. However, if the early disturbance included ions, 
then these tons would not achieve sufficient deflection to miss the 
exit pinhole and any signal attenuation would be very small. 

An unforeseen problem arose in this experiment either because of 
the large distance at which the probe samples the plasma or because of 
the pinhole collimator. The plasma sampled by the probe was not 
reproducible. Major features were always identifiable but their 
amplitudes differed by as mich as a factor of three from shot to shot. 
This unresolved problem has greatly complicated the analysis by causing 
a large deviation in the data. | 

Results of this experiment are reproduced in Figure 17. The straight 
line is a least squares fit to the data. Its slope is 0.19 mV/kV 
indicating that there is little dependence of the early disturbance on 
the applied transverse field. These results suggest that the early 
Signal is not composed entirely of electrons but that there are heavier 


particles present. 


36 





100 


90 
=. Oo o 090 . 
> 80 oo o 
E 00 2 o 
70 O 
o a 
S 60. . ° . 
Q. ° 
© 
e 50 
© r è 
2 40 ® 8 e 
É . i 
30 
20 
10 


2 4 6 8 10 12 14 16 18 20 


Electrode Potential (kV) 


Figure 17. Early Disturbance Peak Amplitude Versus Transverse Electric 
Field Strength. 
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C. EFFECT OF A TRANSVERSE MAGNETIC FIELD 

Because of the complications inherent in the production and inter- 
pretation of data taken with a large transverse electric ibis sel) 
transverse magnetic field was applied to the streaming disturbances. 


The experimental arrangement is depicted in Figure 18. 


target 
normal 





Figure 18. Experimental setup for the application of a transverse 
magnetic field. 

The magnetic field was measured with a gaussmeter. The field in 
the center of the magnetic pole area is 950 gauss. Near the surface of 
CrS leS the field strength increases to 1325 gauss. 

Great care was taken to insure that the incoming laser pulse did 
not strike the target’and that the distance between the magnetic poles 
was sufficient to allow plasma expansion without perturbation. In 
order to determine if the magnetic field was inducing currents in the 
probe a number of pictures were taken with the target removed allowing 
the laser pulse to transit the chamber and exit through port number 


three. Figure 19 shows the oscilloscope response in this situation. 
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Figure 19. Double probe response in the presence of a 950 gauss 
magnetic field. The target has been removed and the 
laser pulse transits the chamber and exits through a 
glass cover over port number three. 

These tracings indicate that no probe response is generated by either 

the magnetic field or laser pulse interaction with components of the 

apparatus other than the target. 


Gyrofrequencies for electrons and ions are given by equation (7): 


(7) o = eB/me (radians/sec) 


The field in this case is 950 gauss so that the ion gyrofrequency 
is 3.35 X 10” rad/sec and the electron gyrofrequeney is 1.67 x 1010 
rad/sec. 

The Larmor radii are given by equation (8) 


(8) Ri = VE 06 


From this expression we find that for electrons the Larmor radius is 
9.75 X 1073 cm and for aluminum ions it is 484 cm. From this it is 
obvious that if the early disturbance is composed essentially of 
electrons, the signal should be greatly diminished at distances greater 


than a few gyro radii or 0.5 cm. However, if the early disturbance is 
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essentially a neutral plasma then there should be little change in the 
amplitude. The ions with their mich larger gyro radius will hardly be 
deflected by the field. As soon as sufficient space charge separation 
is created, the ions will drag the electrons out of the magnetic field. 
Other researches [Refs. 1, 6, 25, 28 and 33] have reported that 

plasma fronts tend to be rich in electrons while the main body is 
ienmbcallyerich. In this case, the front will be changed from electron 
rich to ion rich by the magnetic field. If we consider the plasma on 
a microscopic, interparticle, basis and set the Lorentz force. equal to 
zero we can obtain a measure of the space charge separation necessary 


to pull the electrons out of the magnetic field 


(9) e xXx B = E = VB since Fil B 


(10) Vip B = ke/rº which on solving for r gives: 


1/2 
(11) yr = (ke/VB) where k = 1/47e 


If we substitute Vip = 1.1 X 107 cm/sec and B = 950 gauss into equation 
(11) we find that the separation is 0.1 mm. 

Figure 20 duplicates photographs of the oscilloscopic record of the 
double probe response with and without the transverse magnetic field. 


There is no discernible change in the early signal when the magnetic 


field is introduced. 
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Double probe response at (0, 0, 1.25) to the streaming 
plasma with (a) and without (b) a transverse 950 gauss 


magnetic field. 


eure 20. 


It is interesting to observe that the main plasma is delayed by 


about 10 nsec and the early plasma by about one or two nsec. For the 


main plasma traveling at 10! cm/sec this delay represents a change in 


plasma front position of one millimeter. This compares with the 


prediction obtained from equation (11) if the electron rich area 


originally extends 0.9 mm in front of the neutral boundary of the plasma. 
It might be argued that since the direction of this transverse 


magnetic field is horizontal there is a drift velocity due to 


the gravitational force. Depending on the direction of the fields 


this drift could be either back into Che target or away from the target 


toward the probe. The magnitude of this drift is given by: 
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= ENE 


(12) Vg = me g XB = mcg/qB 
qB 





For the values under consideration in this experiment the magnitude 
of Vog is 5.9 X 107% cm/sec which is insignificant. 

Since the magnetic field has little effect on the early disturbance 
this experiment strongly suggests that the early disturbance is not 
purely electronic in composition. If it were, then the magnetic field 
would deflect the electrons into a gyro radius with diameter much 
smaller than the distances at which the probe measurements were made. 
This would result in the loss of the early disturbance and the probe's 
response would drop significantly. Therefore, if the early disturbance 


is particulate in nature, it must contain heavier ions. 


D. SEPARATION OF THE PHOTOELECTRIC RESPONSE 

During the density contour mapping a disturbance of small half width 
(4 nsec) was observed to occur at approximately t = 16 nsec, i.e. near 
the peak of the laser pulse. The time of occurrence was constant for 
all points on the mapping grid indicating a very high propagation speed. 
Several authors [Refs. 16 and 33] have reported strong X-ray and UV 
radiation emanating from the target area of a laser produced plasma. 
Ready [Ref. 40] has ascribed the early disturbance noted on his 
Langmuir probes to reverse photoelectric effect. His measurements were 
made at 30 cm from the target in near vacuum. As will be noted later, 
the early plasma decays rapidly between 0.75 and 2.0 cm. Therefore at 
30 cm it is most probable that all he saw in terms of early probe 
response was the reverse photoelectric response. At distances less 


than one centimeter from the target the photoelectric response is masked 
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by the much larger response due to the early plasma pulse. Figure 21 
illustrates this masking and shows the appearance of the photoelectric 
Signal at 1.7 and 3.0 cm. 


(b) 


20 nsec/.5V per div. 
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Figure 21. Double probe response to a streaming laser produced plasma 
A OS O) (0. Om lay) andmce) av .(0. 0-3-0). 
Asterisks mark the suspected photoelectric response 

in order to verify that this signal was photoelectric in origin, a 

peed quartz (S10,) plate was obtained. The plate is 2.54 cm square and 

and has a thickness of 0.02 cm. The transmissivity of this plate is 

greater than ten percent for wavelengths longer than 0.12 microns. 

Figure 22 depicts the transmissivity as a function of wavelength for a 

one-half inch plate of the same material. The work function for tungsten 

is approximately 4.55 eV. Therefore, the maximum wavelength for single 

photon photoelectric emission from the probe tips is 0.267 microns. This 


wavelength is marked on Figure 22. The hashed area denotes the trans- 


mitted wavelengths capable of producing the observed photoionization. 
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Figure 22. Transmissivity of a 1/2 inch fused quartz plate. The 
hashed area represents those wavelengths that are 
transmitted and are capable of ionizing tungsten. 

In this experiment, the quartz plate was placed between the probe 
and the target. Probe response with and without the quartz was recorded 


and a comarison made. These oscilloscopic traces are reproduced in 


Figure 23. 





20 nsec / .5V per div. 


— t — 


Figure 23. Oscilloscopic record of a double probe response at (OON 
1.75) with (a) the laser beam passing through the chamber 
but no target in place, (b) without the quartz plate 
screen and (c) with the quartz plate screen in place 
between the target and the probe. An asterisk marks the 
verified photoelectric response 
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This experiment was repeated ten times at different locations in 
the chamber with identical results. As sem,the quartz plate stops 
all physical particles but passes UV radiation which produces the 
Photoelectric response at the probe. 

The height of the response is approximately 300 mV. Assuming that 
en percent Of the electrons ejected fromthe negative electrode are 
collected at the positive electrode this corresponds to a current flow 
of 0.079 amps for 4 nsec or a total of 1.97 X 10? electrons. Assuming 
an average photon energy of six eV and also assuming that ten percent 
of the incident photons result in the ejection of an electron from the 
negative electrode we have an energy deposition on the probe tip of 
1.9 X 1078 joules. The cross sectional area of this tip is 
4.6 X 1073 em. Assuming that the electromagnetic radiation is trans- 
mitted equally in all directions through a hemisphere whose center is 
the focal spot on the plane target, this corresponds to a total energy 
ONA DA 107? joules transmitted at wavelengths between 0.12 and 0.27 
microns. This represents about 1.3 X 10 of the incident energy in the 
laser pulse. 

Equation (25) predicts that for a static plasma the ion density 


mier. 45]: 
(25) a 167 E (m, ART) eh, hema) 


If the 300 mV probe response is due to photoionization of the background 
this would correspond to a background gas ion density of 

IO 7014 em? assuming kT, = 3 eV, m = 28.96 AMU and that the total 
surface area of the probe electrode is 1.4 X 107 ene Since the 
background gas density is only 8.85 X 1011 em? this response must be 


due to the reverse photoelectric effect described by Ready. 
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E. DENSITY MAPPING OF THE EARLY PLASMA 
The expanding plasma density has been mapped using a double probe. 


The mapping grid is shown in Figure 24. 
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Figure 24. Density Mapping Grid 


The general method used was to acquire an oscilloscopic record On 
the probe response at each grid point. A zero time was marked on each 
photograph and corrected for the 35 nsec cable delay. Probe response 
voltages were then read from the photograph corresponding to the 
desired time intervals for the contour maps. These voltages are readable 
to within + 0.2 volts. The voltage from each grid point corresponding 
to a specific time is then plotted on the grid depicted in Figure 24, 
Contours were drawn by hand. All of the data and contour mappings were 
analyzed manually. However, because these are only relative density 
contour maps and because of the many data points involved in each map, 


any errors introduced should be minimal. 
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The area excluded from each map is the corridor through which the 
laser pulse travels in reaching the target. Since plasma density is 
proportional to probe response voltage, these contour maps give 
relative density contours. 

For the early plasma pulse, the probe response was mapped at 10 nsec, 
20 nsec, 30 nsec, 40 nsec and then at 20 nsec intervals to 160 nsec. 
This allows a detailed analysis of the expanding plasma pulse. These 
contour maps are included as Figures 26 to 35. 

The maps indicate that the early disturbance expands rapidly 
(10º cm/sec) along a line inclined at approximately 30 degrees with the 
target normal. This line also describes the path followed by the 
reflected laser pulse. 

It should be noted that typical densities measured throughout these 


maps are on the order of 101° cm 3. 


This is nine orders of magnitude 
below the critical density so that the plasma being examined here is 
essentially transparent to the incoming laser radiation. 

The focal spot size of this laser pulse is 1.5 mm in diameter. It 
should be noted that the early disturbance expands very rapidly to a 
full radius of approximately 1.0 cm by z = 0.2 cm. At z = 0.2 cm the 
plasma radial expansion stops and all further expansion occurs along 
the axis of symmetry. This behavior is strongly suggestive of a force y% 
acting on the plasma to pinch the radial development. Since no 
external fields other than the small probe potential were applied, it 
is evident that this force mst be produced internally within the plasma. 


The self-induced magnetic fields mapped by McKee and Bird [Refs. 8 and 38] 


would provide a pinching effect on the plasma. The J XB force was 
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neglected by McKee as too small to be significant. However, further 
investigation of this force seems warranted and is included in Section 
IV. B- 

It should also be noted that the plasma density gradients are 
higher in the radial direction than in the direction of pulse 
propagation. This is also strongly suggestive of a radial confining 
force. These gradients are plotted in Figure 25. 

Maximum densities noted in the early pulse are approximately 
4 x 1072 cm 73. This corresponds to the 2.5 volt contour in each of 
Figures 26 to 29. In Figure 30 the early plasma pulse has separated 
from the main plasma pulse. Subsequent figures contain only main plasma 
contours. The maximum early plasma densities are greater than the 
background gas density by a factor of five. 

NE monto be noted that as early as t = 10 nsec there is evidence 
of plasma at distances as great as 0.8 cm. This plasma may partially 
be the result of intense background gas photoionization. However, the 
magnitude of the response indicates ion densities three orders of 
magnitude higher than the background gas density (assuming a 
macroscopically quiescent plasma). 

These topographical density maps also show that after flooding into 
tine selaamber , this early disturbance. decays rapidly so that at 100 nsec 
there is little remaining evidence of the early plasma pulse. Another 
interesting feature is that the decay appears to proceed more slowly 
along the target normal than along the axis of propagation of the early 


disturbance (see Figures 30 and 31). 
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Figure 25. Early plasma pulse density gradients in (a) the direction 
of propagation and (b) in the radial direction 
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Figure 26. Laser-produced plasma relative density contours at 
t = 10 nsec. 
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Figure 27. Laser-produced plasma relative density contours at 
w = 2alımsee. 
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Figure 28. Laser-produced plasma relative density contours at 
t = 30 nsec. 
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Figure 29. Laser-produced plasma relative density contours at 
t = 40 nsec. 
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Figure 30. Laser-produced plasma relative density contours at 
t = 60 nsec. 
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Figure 31. 
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Laser-produced plasma relative density contours at 
t = 80 nsec. 
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Figure 32. Laser-produced plasma relative density contours at 
t = J00 nsec. 
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Figure 33. Laser-produced plasma relative density contours at 
t = 120 nsec. 
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Figure 34. 
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Laser-produced plasma relative density contours at 
t = 140 nsec. 
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Figure 35. 
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Laser-produced plasma relative density contours at 
t = 160 nsec. 
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F. DENSITY MAPPING OF THE MAIN PLASMA 

AS an extension of section IV. E., relative density contour maps 
were completed at 100 nsec intervals from 100 to 800 nsec. The same 
mapping grid as shown in Figure 24 was used to produce these maps. 
However, the oscilloscope scales were increased to 100 nsec amd 2 volts 
per division. The results are illustrated in Figures 37 to 44. 

The first, very obvious feature of the main plasma is noted at 60 
nsec in Figure 36. The main plasina body with densities corresponding 
to probe response voltages greater than 14 volts has been shaded. This 
map indicates that at early times Ele main plasma does not expand 
symmetrically about the target normal. The radial density is strongly 
skewed to the left between z = 0.2 and z = 0.8 cm. The high density 
area appears to be centered at (0.6, 7/2, 0.5). Since these contour 
maps were made for only a single plane, no determination of the three 
dimensional distribution is possible. Examination of Figures 30 to 34 
in section IV. E. indicates that this asymmetry decays rapidly so that 
by 120 nsec the main plasma has become essentially symmetric about the 
target normal. To further investigate and hopefully verify the late 
time symmetry of the main plasma with respect to the target normal, a 
special target of 1/8 inch thick 6061 aluminum alloy was constructed 
and mounted onto the face of the half inch target normally used. A 
28 cm converging lens was used to focus the laser pulse to the smallest 
obtainable focal spot size at the target face. This was accomplished 
with an autocollimator. 

Two double probes (numbers one and two) were positioned inside the 


chamber as depicted in Figure 45. 
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Figure 36. Early time asymmetry of the main plasma with respect to 
the target normal 
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Figure 37. Laser-produced plasma relative density contours at 
t = 100 nsec. 
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Figure 38. Laser-produced plasma relative density contours at 


t = 200 nsec. 


63 





RADIAL EXTENT (CM) 





Es 
| ã 
| 2 
300 NSEC 
4 
Figure 39. Laser-produced plasma relative density contours at 


t = 300 nsec. 
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Figure 40. Laser-produced plasma relative density contours at 
t = 400 nsec. 
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Figure 41. Laser-produced plasma relative density contours at 
t = 500 nsec. 
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Figure 42. Laser-produced plasma relative density contours at 
t = 600 nsec. 
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Figure 43. Laser-produced plasma relative density contours at 


t = 700 nsec. 
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Figure 44. Laser-produced plasma relative density contours at 
t = 800 nsec. 
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Figure 45. Experimental arrangement for verification of late time 
plasma expansion symmetry with respect to the target 
normal. Probe one is at (0.2, 7/2, 2.5) and probe two 
A Sa oe 

A total of 50 shots were made without disturbing the target. 

Pictures of the double probe responses were taken every five shots as 

displayed on a Tektronix 555 dual beam oscilloscope. The laser had 

burned a visible hole into the target after the fourth shot. 

Reproductions of the oscilloscopic traces at ten shot intervals are 

included as Figure 46. 

Some increase in the signal at probe number two was expected due 

to a shotgun effect as the plasma expands from the base of the deepening 

hole. However, as the traces in Figure 46 indicate, there is no 

appreciable change in the plasma density sampled by either probe during 
the entire sequence except that the early plasma pulse is somewhat 
attenuated after the fifth shot. This behavior strongly suggests that 


the expansion dynamics are independent of the plasma behavior prior to 


its rapid expansion into the background. 
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Figure 46. Series of oscilloscopic traces taken on a dual beam 
oscilloscope during the burning of a hole 1.029 cm deep 
in an aluminum target. (a) probe two and (b) probe one 


t 


Following this experiment, a detailed examination of the crater was 
undertaken with an electron microscope. Figure 47 is an electron 
micrograph of the entry point of the laser pulse, Figure 48 depicts 
the exit hole from the 1/8 inch plate, Figure 49 shows a cross section 
of the hole in the 1/8 inch plate, Figure 50 is a cross section of the 
crater burned into the 1/2 inch back plate, Figure 51 shows a cross 
section of the base of the hole excavated into the back plate and 
Figure S2a@eer1cts a typical crater wall. 

Careful measurement of the volume of the cavity reveals that 
approximately 0-014 cm of aluminum were excavated. This represents a 
total of 8.44 x 10% aluminum atoms or approximately 1.7 X 1019 atoms 
per laser pulse. This value is two orders of magnitude larger than 
the 1017 atoms reported by Honig [Ref. 23]. However, examination of 


Figure 47 indicates evidence of splashes of mojten material radiating 
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Figure 49. Cross section of the hol® durru'd through the 1/8 inch 


aluminum target plate. ragnification is 231. 





Figure 50. Cross section of the crater burned into the 1/2 inch 
aluminum back plate. Magnification is 17X. 
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BOLOS 51. Base of the crater burned into the 1/2 inch aluminum tach 
plate. Magnification is 480X. 





Figure 52. Typical section of the crater wall in the 1/2 inch aluminum 
back plate. Magnification is 540X. 





from the entry point of the laser pulse. Examination of the probes 
after the experiment reveals that each is heavily coated with a layer 
of shiny aluminum. These splashes and the aluminum coating have not 
been observed in previous experiments using the laser in a defocused 
mode. This evidence indicates that a great deal more material was 
ejected during the hole burning process than is ejected when the laser 
pulse is defocused to produce only a shallow crater. Rather than being 
ionized, this additional material appears to be essentially neutral and 
may even be in a partially condensed state. This statement is 
supported by lib coating on probe number two which did not 
indicate the presence of a significant number of ions at any time during 
the experiment. 

An obvious cause for the increase in ejected material is that the 
hole tends to contain all of the incident laser energy whereas in the 
shallow crater produced by the defocused laser, a significant amount 
of the incident energy may be reflected from the target and expanding 
plasma. Examination of Figure 51 reveals evidence of dynamic shock 
processes in the base of the crater. large pieces of aluminum are 
apparently being ejected by shock waves. The ragged appearance of the 
entire surface of the crater suggests that dynamic removal of material 
may be more important than thermal processes. Inhomogeneities in the 
aluminum alloy structure may strongly contribute to the appearance. 
This subject will not receive further consideration in this 
investigation. 

Main plasma density gradients in the radial and axial direction are 
plotted in Figure 53 for times of 100 nsec, 400 nsec and 800 nsec. The 
radial density gradient is plotted at z=1.0 and z= 2.0 cm. The axial 


density gradients are along the target normal. 


ie 





These main plasma density contours lend further evidence for radial 
confinement of the plasma. This evidence is seen in the cylindrical shape 
of the expanding plasma out to 3.0 cm (see Figure 40) and by the 
initially steep and slowly diffusing radial density gradients seen in 
Figure 53. 

The high density gradients for small z values at long time intervals 
seen in Figures 43 and 44 must be the result of rapid plasma cooling in 
contact with the target face since the laser pulse shuts off at 
t = 30 nsec. 

To obtain a more detailed picture of the plasma density gradients 
close to the target under slightly different conditions, data was collected 
on a one millimeter grid from z = 0.2 cm to z = 1.2 cm and radially to 
1.5 cm. Contour maps were constructed at 100 nsec intervals from 200 nsec 
to 800 nsec. These maps are included as Figures 55 and 56. The data were 
obtained with probe number three biased at -15 volts and with the laser 
Ensine reduced to 4.0 + 0.4 J. 

The most striking feature of the plasma is that it retains nearly 
the same density gradients and radial dimension throughout the time scale 
of this survey. This is most dramatically illustrated in Figure 54 where 
the one volt contour for 200 nsec is superimposed on the one volt contour 
for 800 nsec. As can be seen in Figure 54, the radial extent has 
increased from 0.7 cm to only 1.0 cm during this time interval. This 
provides further evidence for the presence of a radially confining force. 
Also note that between 0.5 cm and 1.2 cm the density gradients are almost 
entirely in the negative radial direction so that any temperature gradient 
in the negative axial direction would provide the necessary source term 


for the reversed self-induced magnetic field reported by Bird. 
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Figure 53. Main plasma density gradients at (a) 100 nsec, (b) 400 nsec 
and (c) 800 nsec. 
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Figure 54. Density contours for (a) 200 nsec and (b) for 800 nsec. 
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Laser produced plasma density contours at (a) 200 nsec, (b) 
300 nsec, (c) 400 nsec, (d) 500 nsec and (e) 600 nsec. 
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Figure 56. Laser produced plasma density contours at (a) 700 nsec 
and (b) 800 nsec. 
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G. DEPENDENCE OF THE EARLY PLASMA ON LASER PULSE ENERGY 

The maximm density in the early sen has been found to vary 
linearly with the laser pulse energy. A series of oscilloscopic traces 
of the double probe response at (0, 0, 1.9) were recorded with the 
laser pulse energy varying from 2.9 joules to 9.2 joules. The peak 


early probe response versus laser energy is plotted in Figure 57. 


ho 


Early Peak Probe Response (v) 





6 7 
Laser Energy (Joules) 


Peak Double Probe Response to the Early Plasma Versus Laser 
Pulse Energy 


It appears that the early response has a linear dependence on the 


Figure S. 


laser energy. For this particular probe (probe number three) and 


within the energy range plotted, the dependence can be expressed in the 


relation: 


(13) Vop = (1.18, - 0.375) volts 


Where Vpp is the peak probe response to the early plasma pulse and Er, 


is the laser pulse energy in joules. It should again be noted that these 
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data and equation (13) are valid only for this experiment. The 
amplitude of the early pulse is highly dependent on background gas 
pressure as shown by Arifov [Ref. 5]. These data were collected in 
near vacuum at 107? torr. This pressure is at the threshold of the 
pressure dependent region described by Arifov. There are peak ion 
densities of 2.5 X 10!º em 73 in the early plasma pulse at this pressure. 
As the early plasma pulse propagates, it ionizes the background. 
However, if the background density is less than the ejected plasma 
density then the additional ionization will be comparatively small. 
The background density is equal to the early plasma density ve 
pressure of 7 X 107? torr. As the background pressure increases above 
this value, ionization by the streaming plasma increases the total ion 
density detected by the probe and the pressure dependent region is 


> 


reached. The significant boundary at 7 X 10“ torr is consistent with 


Arifov's and Isenor's work [Refs. 5 and 25]. 


H. TIMING SEQUENCE FOR DISTURBANCES LEAVING THE TARGET 

An attempt has been made to determine the times at which the various 
disturbances left the target. This was accomplished by assuming that 
the speed of the disturbance is constant between the target and the data 
point. By noting the time at which the disturbance reaches a specific 
data point and subtracting the transit time using the measured speed of 
the disturbance an estimate of the time at which the disturbance left 
the target is obtained. This method appears to give reasonable results 
for the early disturbance and the photoelectric response due to the 
constancy of their propagation speeds. Repeated measurements with 
available data were averaged. The average computed departure time is 


t = 16 + 2 nsec for the photoelectric disturbance and t = 18 + 6 nsec 
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for the early plasma disturbance. The average of the times at which 
the main plasma left the target is t = 2 + 15 nsec assuming a constant 
speed of 1.1 X 10! em/sec. This last value does not seem consistent 
with any dynamic model for the plasma generation. It is assumed that 
the error is introduced because the main plasma speeds are much higher 
close to the target. It is not possible to measure the position of 
the main plasma front at distances less than 0.6 cm because the early 
disturbance is still separating from the main plasma body. This is 
illustrated in Figure 58. This problem plus the higher oscilloscope 
time scale (100 nsec per division) used to measure the main plasma 


account for the high standard deviation in this computation. 









Main Plasma 


probe response 
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—Time——> 
Early Disturbance 





Figure 58. Double probe response at (0, 0, 0.6) showing incomplete 
separation of the early disturbance from the main plasma 


body. 
I. CORRELATION OF THE QUARTZ PRESSURE PROBE AND DOUBLE PROBE RESPONSES 
In this experiment a double probe and a quartz pressure probe were 
placed side by side at an axial distance of 1.4 cm from the target. 
Probe responses were recorded from each diagnostic probe on successive 


shots using the same delay cable. 
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In this investigation, the momentum transport of the plasma provides 
the strain on the quartz pressure probe. Assuming that plasma particle 
collisions with the pressure probe are elastic and that the mass of the 
plasma is much less (6 orders of magnitude at least) than that of the 
probe and its holder, the force density producing the strain in the 
gauge is: 


(13) F= = E shy + mvn, + mv an) 


Since the electron mass is negligible compared with the atomic mass we 


have m = m > m, and equation (13) can be approximated by: 


(14) Fs sen TP (n, + n,)) 


It should be noted that this equation assumes that the plasma flow 
velocity, VI p>» is much larger than the plasma ion thermal speed. This 
assumption is justified in Appendix A. 

Applying the assumption of elastic collisions and doing a simple 
dimensional anaylsis on equation (14) indicates that: 


(15) F= AP = OmV.,(n. + AL 
= BR Dune 


Where Tp is the time for the plasma to travel one cm. Since 


= (V. 2 equation (15) can be rewritten as: 


(16) F = 2m(V. mae (n; + N) re 


This analysis describes the plasma as a square pulse which it is 
not. However, for the purposes of this discussion the approximation 
of a square plasma pulse will do. For the main sen pulse with 
m = 27 AMU, Vip = 1.1 X 107 cm/sec and (n; +n,) = 2 X 1019 om, the 


strain producing force density is 0.08 newtons/cm3. 
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A similar analysis of the early plasma pulse with Vip =1.1X 108 
cm/sec and (ni + nn) WEDER 1012 cm” yields a strain producing force 
Gensity of 0.01 newtons/cm. 

The ratio of these two stress producing force densities is 0.125. 
This is also the ratio of voltage outputs from the crystal in the 
presence of each pulse since the probe response is proportional to the 
stress. | 

Ihese impulsive forces produce stress waves in the aluminum gauge. 
These waves propagate with the speed of sound. For aluminum, 

U, = 6,420 m/sec [Ref. 4]. Therefore it takes 390 nsec for the stress 
wave to propagate through the 0.25 cm aluminum cover disc and epoxy to 
the quartz gauge. 

Figure 59 depicts the signals produced from the same position by 
(a) the double probe and (b) the pressure probe. As can be seen, the 
main plasma front corresponds very well with the large, late response 
from the pressure probe. However, the early response from the pressure 
probe occurs approximately 200 nsec before a stress wave generated by 
the early plasma would arrive at the crystal. No conclusive explanation 
is offered for this temporal discrepancy in the arrival of the early 
pressure probe response. Perhaps the speed and sharpness of the early 
plasma cause the wave to travel through the aluminum at a speed in 
Bess of the speed of sound. The electrodes are well shielded with 
epoxy and the outer case is isolated from the circuit so that electro- 
static pickup is unlikely. However, the early signal does correlate 
with the time an expected electrostatic pickup would be measured. It 


is also possible that the obvious early signal is not that associated 
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with the early momentum transfer. The expected amplitude of the early 
signal is approximately 4 mV. A signal of this magnitude could easily 
be lost in the noise. 

The fact that there are two distinct disturbances registered by the 
pressure probe indicates the presence of two momentum transporting pulses. 
However, the poor temporal correlation of the early signal makes the 
results of this experiment somewhat inconclusive. É 

It should be noted that researchers at the Stanford Research 
Institute have fired high power, 30 nsec, laser pulses at the gauge's 
aluminum target face. The quartz gauge then measures the momentum 
transported to the target by the fast blow off of laser produced plasma. 
These investigations have revealed two momentum impulses delivered to 
the target as depicted in Figure 60. 

(a) 


100 nsec/500 mV per div. 


double probe 
response 


pressure probe response 
| 4 | target removed 


b) 100 nsec/20 mV per div. 


a | 
PITA 
em main plasma impulse expected 
time early impulse expected 


pressure probe 
response 


early pressure probe response 


Figure 59. Correlation of signals received by (a) the double probe 
and (b) the pressure probe at (0, 0, 1.4) in the presence 
of a laser plasma. 
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Figure 60. Stress time history in a laser irradiated sample in 
vacuum. r* is the propagation time of the stress 
pulse through the aluminum cover. r is the 
propagation time of the stress pulse through the 


quartz guage. 

J. AZIMUTHAL MAGNETIC FIELD SYMVETRY 

To verify the azimuthal symmetry of the self-induced magnetic field 
a series of pictures were taken at 9 = 5° intervals at (O Cno). 
The magnetic probe was inserted through a holder in port five of the 
chamber (See Figure 4). The probe holder is circular and nests in a 
circular flange allowing rotation of the probe without disturbing the 
r or Z coordinates. 

Figure 61 reproduces the probe response at 15 degree intervals from 
6="7 to 8 = 7/2. The gross aspects of the magnetic probe response 


are symmetric. 
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Figure 61. Magnetic Field Symmetry for (1.0, 9, 1.5). Scales are 
100 nsec and 100 mV per division. 


—Time—> 


Unfortunately this experiment was conducted prior to discovery of 
the plasma expansion asymmetries. Therefore the oscilloscope scales 
were 100 nsec and 100 mV per cm. These scales and the z coordinate of 
1.5 em do not allow resolution of the early magnetic field asymmetries 
expected with the plasma density anisotropy. Therefore, this 
experiment should be repeated with r = 0.5 cm, 9 = zero to nm and 
z = 0.5 cm. The oscilloscope scales should be set at 100 mV and 


Wenn ee per division. 
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IV. ANALYSIS OF EXPERIMENTAL RESULTS 


In this section two facets of the experimental results are examined. 
In section (A), a comparison is made between the particle content in 
the early plasma pulse and the number of atoms initially exposed to the 
laser radiation (to a depth equal to the skin depth). In section (B) 
the J XB force density is analyzed as a possible source for the radial 
confinement. 


A. ORIGIN OF THE EARLY PLASMA PULSE AND ITS PRESSURE DEPENDENCE 
Assuming a macroscopically neutral plasma and plasma symmetry about 
the axis of propagation, the total number of ions involved in the early 
plasma pulse is approximately 2 X 1013 Al ions. The skin depth for 
aluminum is 61 angstroms. The laser focal spot size is 0.15 cm. The 
total volume initially exposed to the laser radiation is 1.08 X 1078 cm’. 
This corresponds to a total of 6.5 X 10!" aluminum atoms. This agrees 
to within one order of magnitude with the number of observed ions. If 
the early plasma pulse also has neutrals by the time it reaches the 
probe grid then this agreement will be even better. Therefore, at least 
in terms of total numbers, the early blowoff of plasma could be due to 
those atoms which are directly stimulated by the initial laser radiation. 
The pressure dependence of the early plasma has been previously 
described by Arifov [Ref. 5]. Assuming that the early blow off of 
material is independent of background gas density implies that there 
are initially 2.4 X 1012 ions/cm in the fast stream. Background pressures 
reach this density at 7 X 10 torr. This pressure is consistent with 


onset of the pressure dependent region. The fast streaming early plasma 
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further ionizes this background. However, when the background pressure 
is low, the additional ions do not add appreciably to the total. As 
background densities increase above the threshold there are 
Significantly more ions created by the fast streaming plasma. As the 
background pressure is increased above 250 mtorr the fast streaming 
early plasma mist be stopped in a short distance causing a reduction 


in early signal strength. 


B. ANALYSIS OF THE J XB FORCE DENSITY 

McKee [Ref. 38] using typical values for plasma expansion into a 
250 mtorr background found that the J XB force density was much less 
than the radial pressure gradient force density. McKee also plotted 
this force density at 120 nsec for 250 mtorr of nitrogen background gas. 
He found that the magnitude of the largest force at this time is 
et x 107 newtons/ 3, À reproduetion of his computer produced plot is 
included as Figure 62. 

Bird [Ref. 8] has found magnetic fields of the order of 200 gauss 
at (0.4, 1/2, 0.75) for aluminum plasmas. This measurement was made at 
200 nsec after arrival of the laser pulse. The pressure was 0.1 microns 
which is close to the pressure independent region for both the magnetic 
field magnitude and magnitude of the early disturbance as reported by 
Arifov [Ref. 5]. The current corresponding to this 200 gauss field is 
of the order of 400 amps/cm*. Using maximum current and magnetic field 
densities indicates the presence of a J XB force density of approximately 
7 X 10% nt/m. This is consistent with McKee's maximum values. 

At 500 nsec, the extent of the plasma (2 volt contour) is 


approximately r = 1.5 cm by z = 4.5 cm. Assuming symmetry with respect 
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Figure 62. J XB force an at 120 nsec for 250 mtorr of nitrogen 
background gas. SER a of the largest force at 
this time is 1.4 X 107 nt/m3 
to the target normal this corresponds to a total volume of 32 cm. If 
1017 atoms are ejected and ten percent are ionized, then the ion density 
would be on the order of 3 X 1014 ions/em. 
In section ILI. F. ion densities for the streaming min plasma of 
NR x 1045 Vp (Vp is the probe response voltage) were found. Typical 
probe response voltages are on the order of ten volts for the main plasm. 
This corresponds to a total ion density of 1.4 X 1014 ions per om, 
These two results are in good agreement. The total plasm density 


including neutrals is then taken to be 1012 cm3. 
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Assuming a plasma temperature of three electron volts, and a 
characteristic density gradient length of f = 0.8 cm, the pressure 


gradient force density becomes: 
(17) nkT/P = 10° newtons/m3 


| Using either the maximum J X B force density found by McKee or that 
found on page 90 we see that the J x B force density is of the same 
order of magnitude as the pressure force density found in equation (17). 
Therefore the J XB force density is capable of at least partial 
confinement of the plasm and in any case should be taken into account 
in describing the dynamics of the expanding laser produced plasm. 

in order for this confinement to continue for the ennen times, 
the characteristic diffusion times must be at least of the same order 
as the observed effect. Otherwise the plasma will diffuse through the 
magnetic field. 


The magnetic field diffusion time is given by: 
i 


(18) + = pof 


m 


The characteristic length ( f ) of the magnetic field at 200 nsec is 
estimated from McKee's data to be 0.8 cm. The D.C. conductivity of the 
plasma is given by: 


(19) 0 = ne“ /m, Vo 


where uv, is the electron collision frequency. For electron-ion and 
electron-neutron collisions the collision frequencies are given by 


Tanenbaum [Ref. 46]. 
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1 
(20) vi 


3.62 X 107 n,7,2/2mA sec. 


ea) 


Dos 
il 


2.60 X 10, m7 Y ea? see ~l 


Where T is in °K, n in particles/m, Ina is typically about ten, o is 
of the order 101%. assuming T, = 4 eV, n, = 1020 m3 and n, = 102% m3 
the collision frequencies are found to be Vej = 3.62 X 108 sect and 
DE DGA 10! sect. Substituting these values into equation (19) 
reveals a D.C. conductivity of 7.8 X 103 mo/m. Using this value in 
equation (18) we find that the magnetic field diffusion time is of the 


order: 
(22) ee 626 nsec 


This value appears consistent with McKee's work and explains how the 
plasma could be radially confined within the time scale surveyed in this 
investigation (t< 800 nsec). While not plotted, an examination of the 
data indicates that the plasma becomes very diffuse after 900 nsec. 

This is due in part to diffusion but is more probably due to general 


cooling and recombination. 
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V. SUMMARY AND CONCLUSIONS 


The early disturbances have been separated into "reverse" photo- 
electric response due to photoemission of electrons from the probe and 
a series of at least two additional disturbances traveling with constant 
speeds of 1.1 X 108 cm/sec and 5.9 X 10! em/sec. The first of these 
disturbances has been analyzed in detail. 

The early disturbance is found to propagate along a line describing 
the path of the reflected laser pulse. The disturbance is symmetric 
about this propagation axis. 

After separating the photoelectric response, consideration must be 
given to the possibility that the remaining disturbance is created by an 
ionizing potential wave or a blast wave. References 2, 9, 22 and 47 
provide some insight into the ionizing potential wave phenomenon. In 
this experiment, the background gas density is approximately 9 X 1011 em. 
Assuming that the temperature of the background is three electron volts 
after ionization by the potential wave we can use equation (25) from 
Appendix A to obtain an estimate of the quiescent plasma ion density 


necessary to generate the observed probe current of 0.66 amperes at 


OO). 


We 


(26) u 1867 I, (m, /KT,) /eh, (E) 


The average molecular weight of air is m = 28.96 AMU = mj. The surface 
area of the electrodes is 0.0146 em*. Inserting these values into 
equation (25) indicates that the required ion density is 1.5 X 1015 ion 
charges/cm3. This ion density is three orders of magnitude greater than 


the background gas density. This argument appears sufficient to exclude 


94 





lonization of the background gas by any means as a possible source for 
the early disturbance. 

Ihe higher densities observed can only be supplied from the target 
in the form of a fast blowoff of material. This possibility has been 
considered by many researchers [Refs. 10, 15, and 24]. However, a 
review of the literature indicates that this early blowoff is generally 
thought to be composed entirely of electrons. The results of the present 
research in which transverse electric and magnetic fields were introduced 
indicate that the early disturbance is composed of ions and electrons. 
Ion densities of 2.4 x 1014 cm73 Bene Bean measured in this fast streaming 
plasma. Electron densities may be slightly higher than this creating a 
macroscopically negative plasma potential. This would explain the 
previous observations. 

The pressure dependence of the early plasma has been previously 
described by Arifov [Ref. 5]. Assuming that the early blowoff of material 
is independent of background gas density implies that there are initially 
2.4 X 1012 ions/cm? in the fast stream. Background pressures reach this 
density at 7 X 107° torr. This pressure is consistent with onset of the 
pressure dependent region. The fast streaming early plasma further 
lonizes the background. However, when the background pressure is low, 
the additional ions do not add appreciably to the total. As background 
densities increase above the threshold there are significantly more ions 
created by the fast streaming plasma. As the background pressure is 
increased above 250 mtorr, the fast streaming early plasma must be 
stopped in a short distance causing a reduction in early signal strength. 
In this investigation the early plasma density is strongly attenuated 


between 0.8 and 2.0 centimeters. However, a small signal is still 
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observable out to the limits of this survey (5 cm). Residual 
ionization decreases rapidly so that by 100 nsec there is little 
remaining evidence of the passage of this early pulse. 

The main plasma speed has been measured at 1.1 X 10/ cm/sec. Early 
time anisotropies and asymmetries are observed in the expansion of the 
main plasma. This results in an initial high density area on the side 
of the target away from the incident wiser ace This asymmetry is 
short lived and for times greater than 120 nsec, the plasma can be 
described as symmetric about the target normal. This late time 
Symmetry in the free expansion is evident even when the plasma is created 
in a deep hole and gains momentum in a direction off the target normal 
before reaching the surface of the hole. In this investigation a one 
centimeter deep hole was excavated into an aluminum alloy target. 
Throughout the excavation there was little observed change in the plasma 
expansion. This suggests that the expansion dynamics are independent 
of the manner in which the plasma reaches the surface of the target. 

Previous investigators [Ref. 23] have observed approximately 1017 
particles ejected per incident laser pulse under similar conditions. 

In this hole drilling exercise approximately 1019 atoms were ejected per 
incident laser pulse. Most of this material was neutral and completely 
coated the exposed probes. The reason for identifying the additional 
material as neutral is that little change was observed in the ion 
densities at probe one and probe two though heavily coated with aluminum, 
never did register a significant ion density. 

Examination of the cross sectioned crater indicates that a 


significant amount of the target damage is generated by dynamic shock 
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waves rather than by thermal processes. This is especially true of the 
base of the crater where whole chunks of material are being ejected, 
possibly by dynamic shock processes. 

The main laser plasma is observed to expand very rapidly in the 
radial direction to a diameter of approximately 2.4 cm. This maximum 
radial extent is achieved at an axial distance of 053 emo At this Beine 
the plasma plume ceases to grow in the radial direction and expands 
symmetrically along the target normal with a cylindrical shape. Peak 
plasma densities in the main plasm pulse are approximately 
2.5 X 10l” cm3. Radial density gradients remain high for at least one 
microsecond while the axial density gradient decreases steadily with 
increasing time. This behavior is strongly suggestive of the presence 
of a radial confining force acting on the expanding plasma. The self- 
induced magnetic fields generate a J XB force density which is seen to 
be comparable (within an order of magnitude) with the radial pressure 
gradient force density. This confining force is postulated as an 
important element in creating the observed cylindrical shape of the 
laser produced plasma. This D X El force must be included in any dynamic 
model of the expanding plasma. 

Table (3) summarizes the properties of the various plasma features. 


PLASMA FEATURE SPEED ION CHARGE TOTAL TOTAL ENERGY 
(cm/sec) ¡DENSITY (cm”3), NUMBERS (joules) 


— 320x100 - | | 8x208 
First Barly Plasma 0.55 
Second Early Plasma 1 x 1013} 0.079 
Main Plasma 1.1 
AAA 
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The total energy radiated by photons with wavelengths between 
0.12u and 0.267u has been previously computed. The energy content in 
each of the two early plasmas assumes that each of the plasma particles 
has a directed speed equal to the speed of the plasma front. The 
ionization energy and assumed early plasma temperature of ten electron 
volts contribute a negligible amount of energy. 

The energy of the main plasma is much more difficult to compute. 
Using the methods outlined by Dawson [Ref. 14] the temperature in the 
main plasma body is estimated to be T = M, V“ p/5(Z + 1). Early in the 
expansion, the average ionic charge is assumed to be Z = 3. This gives 
a plasma temperature of 170 eV and a total plasma energy of 1.1 joules. 

Examination of Figures 37 to 44 indicates that only a small portion 
of the particles in the main plasma have the directed velocity seen in 
the plasma front. The bulk of the main plasma does not move from the 
immediate vicinity of the target. The sum of those energies accounted 
for in Table (3) is 3.2 joules or about 65 percent of the energy in the 
incident laser beam. The energy absorbed in the target by conduction 
and shock wave absorption has been estimated at 20 to 30 percent of the 
incident laser pulse energy by Peter Krehl at the Stanford Research 
Institute (personal communication). 

This account of the incident energy neglects energy contained in the 
neutral particles which constitute 90 percent of the ejected material, 
reflected laser pulse energy and plasma radiation at frequencies other 
than those between 0.12 and 0.267 microns. These mechanisms must 


account for the additional 35 percent of the incident energy. 
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VI. AREAS WARRANTING FURTHER RESEARCH 


Because of the early time asymmetry of the expanding laser plasma, 

a three-dimensional map of both the early and main pulses is necessary 
to fully understand their extent. In conjunction with this, the data 
obtained in this investigation should be used to provide a computer 
generated map of the plasma density. 

Development of plasma temperature diagnostic techniques at NPS would 
allow production of similar maps for the plasma's isotherms. With this 
information, the school's computer should be used to plot the cross 
product of density and temperature gradients. This would predict the 


magnitude and direction of the magnetic field due to the VkT, X Vn, 


= Ne 








source term described by McKee and Bird. These results could then be 
compared with McKee's measurements [Ref. 38] in the pressure independent 
region described by Bird [Ref. 8]. 

This investigation has revealed asymmetries in the early expansion 
of the laser produced plasma from a plane target. The dependence of 
these asymmetries on the angle of incidence of the laser pulse requires 
investigation. In conjunction with this, an accurate measure of the 
reflected laser pulse energy at 1.06 microns should be correlated with 
the angle of incidence. 

The pressure dependence of these early time asymmetries could help 
explain the pressure dependence of the early plasma density and 


magnitude of the self-induced magnetic field strength [Refs. 5 and 8]. 
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Measurement of currents in the target generated by the inhomogeneous 
blowoff of ions and electrons and by the self-induced magnetic fields 
could help resolve the induced magnetic field and plasma density 


characteristics at distances too close to the solid target material to 


be measured with external probes. 
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APPENDIX A 


I. DOUBLE PROBE THEORY 


When a potential difference is applied between two electrodes 
immersed in an ionized but macroscopically neutral plasma, the positive 
ions are attracted to the negative electrode and the plasma's electrons 
are attracted to the positive electrode. This action tends to destroy 
the plasma's > Dy creating space chre Selectric fields: ie 
plasma responds ‘to this perturbation by creating a sheath which restricts 
the charge separation to a short distance from the intruding electrodes. 
This sheath absorbs the potential difference in a very short distance 
which is of the order of a few Debye lengths. The Debye length is equal 
to: 


(23) A = 69X (n/a 37° (meters) 


Assuming a temperature of 100 eV and an electron density of or 


cm”, the Debye length is 7.4 X 107% cm. The bulk of the 
plasma remains charge neutral. 

When the surface area of the two electrodes is equal and the entire 
system (electrodes, biasing power supply, scopes, etc.) is isolated from 
ground it is said to "float" with the plasma potential. 

If the biasing potential is initially zero and the double probe is 
immersed in a neutral plasma with T; = di then the random electron flux 
to the probe tips will be larger than the random ion flux due to the 
higher electronic thermal velocity. This will cause the probe tips to 
collect a net negative charge until the potential change is just 


sufficient to allow an equal flux of ions and electrons. ‘The probe tips 


TOT 





will then be at a potential below that of the plasma. This potential is 


called the floating potential. This is illustrated in Figure 63. 


Probe Potential 


SD ng 





Electrode } Electrode 2 


Figure 63. Probe potential with zero bias applied. 


If the number two probe is now biased slightly negative, Vo will 
become more negative and Y, will become slightly less negative. Thus 
more electrons will flow to probe one. This results in a positive 
current flow to probe two. When the number two probe potential becomes 
negative enough to exclude electrons and collect the saturation ion 
current, probe number one will be at fa potential just sufficient to 


allow collection of the net electron flux necessary to cancel the ion 


current to probe one. This condiction satisfies Kirchoff's Law: 


A) | Tail u A [Lo | 
The potentials involved in this situation are depicted in Figure 64 and 
the current flows in Figure 65. 
An important advantage of the double probe is that the current in 
the system is limited by the ion saturation current. In the absence or 
imbalance of ions, the potential of the entire system changes so that the 


maximum current drawn is always limited by the ion densities present. 
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ion current at probe number two. 


n, = 1.671, (m,/et,)1/2/eA, (en) 
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Collection of saturation ion current at probe two. 


PROBE 72 

pe — 

a he Na+ 
Currents flowing in a double probe collecting saturation 


For a macroscopically static plasma, the saturation ion flux to a 


cylindrical probe of surface area (A) is given by [Ref. 45]: 








The electron temperature co can either be estimated or found from 


the slope of the probe characteristic at the origin [Ref. 11). 


(26) aI/dV = e (14) å (154) 
Le THE 
1+ 2+ 
If the areas of the two electrodes are equal and the product of the 
probe bias voltage and the average ionic charge is much less than the 
average thermal energy of the ions, then 1, el 
into equation (26) and solving for T, we find that: 


24" Substituting this 


(27) T, = el, (1/al/av)/2k 


It should be remembered that these equations ame valid only ina 
Static plasma. In the current investigation, the initial random thermal 
energy is rapidly converted to ordered flow kinetic energy. This is at 
least true for plasma expansion into near vacuum conditions. Dawson 
[Ref. 14] finds that the expansion appears to have taken place by the 
time the plasma has expanded to dimensions of 0.01 to 0.1 em. If the 
flowing plasma is to remain charge neutral, then n, = ns Zs and the flow 
velocities of ions and electrons mst be equal. Therefore, the ordered 
kinetic energy of the plasma is ar due t20 the ions because 
ASA greater mass. 

Assuming that T, = T, = 100 eV and a flow velocity of 107 Cm SEC, 
the ratio of ion thermal to ion directed speed is 2.5 X 107 and the 
ratio of the electron thermal speed to flow speed is approximately 13. 
Therefore the smaller current will be that of the ions. 

Now I, = V/R where R is the resistance across which the current is 
measured. This is also the average number of positive charges collected 


by the probe per unit time. This can be expressed as: 
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(28) WR= MAND 


Where n;Z, is the average electronic charge per unit volume and (A) is 
the cross sectional area of the probe perpendicular to the plasma stream 
velocity (Vip): 

Equation (28) does not include secondary emission at the probe tips 
due to ion bombardment. Koopman [Ref. 30] has measured this effect and 
finds that the secondary emission coefficient varies between zero and 
0.5 electrons emitted per incident ion. He also found that secondary 
emission was not strongly dependent on the ion kinetic energy in the 
range of energies encountered here. Therefore, the secondary emission 
can be treated as a constant. To include this effect, equation (28) 
can be rewritten to include an average number of secondary emissions 


per incident ion: 


f NEDRE 


Me 
(29) V/R = n,2, (1 +e) A Vip 


Since A, V, R and Vip are measured quantities and e is a constant, 
equation (29) provides a reasonable estimate of the absolute plasma 
density. For equation (29) to be valid we require that the probe 
electrodes be of equal area and that the plasma be homogeneous across 
the probe tips. The probes were carefully constructed to satisfy the 
first condition. Previous workers at NPS have been careful to keep 

the plane containing the probe tips perpendicular to the direction of 
the streaming plasma. However, the probe tips are about one millimeter 
apart and the characteristic density gradients have lengths of five 
millimeters. Therefore, it was suspected that this orientation was not 


critical. To test this hypothesis, the probe was maintained in the 
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same geographical position and rotated through a variety of angles. 
Figure 60 shows the results as recorded from the oscilloscopic tracings. 
These reproductions show that the signals are essentially unaffected by 
rotation of the plane containing the probe tips. 

Qn the basis of this experiment, the plane containing the probe 
tips was maintained perpendicular to the z-axis throughout the 


experiments regardless of the expected plasma flow velocity. 


— $ — 


(a) (b) 


Figure 66. Effect of rotating the plane containing the probe tips with 
respect to the targets z-axis. The angles (Ø) given are 
the angles between the plane defining the probe tips and 
the target normal. (a) recorded at two volts and 100 nsec 
per division and (b) is recorded at 500 mV and 20 nsec per 
division. 
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IT. DOUBLE PROBE CONSTRUCTION 


Figure 67 depicts the probe circuit used throughout this investigation. 
In one instance the four microfarad capacitor was replaced with a 0.5 
microfarad capacitor. No detectable difference in the probe response was 


noted. 


Electrodes 





stainless steel 


y shield ground 


Power | 
“RE -BNC A O- O 


GL eo rm J 


resistor 1079 
$- A lo 


aluminum housing Capacitor 


Figure 67. Double Probe Circuit 


107 





The following tests have proven effective to insure proper probe 


operation. 


1. Check each electrode for continuity (no resistance between 
(A) and (I) and between (B) and (J)). 


2. Insure isolation of the probe circuit from the aluminum housing. 
(>100 kQ between (H) and (E) and between (C) and (E)). 


3. Insure adequate grounding of the stainless steel probe shield. 
(zero resistance (F) to (E)). 


4. Check the current limiting resistor. (>5k2(G) to (B)) 
>. Check the current viewing resistor for one ohm. ((D) to (C)). 


6. Check the circuit for continuity. (No resistance (H) to (D) 
and one ohm from (H) to (C)). 


Most of the glass enclosed probes were made by a professional glass 
blower. These probes had very fine tungsten electrodes with a diameter 
of 0.013 cm. The electrodes were 0.25 cm long and were set approximately 
0.08 cm apart. The diameter of the glass tip was approximately 0.1 cm. 

At one point in the investigation, a larger signal to noise ratio 
was desired. This was necessary when examining the early disturbance at 
Gistances greater than two centimeters from the target. In order to 
achieve this higher ratio a double probe was constructed with large 
electrodes. Production of this probe required little technical expertise 
and the method is repeated here as it requires a minimm of equipment and 
could be duplicated in any laboratory. 

The electrodes were constructed from 0.05 cm diameter tungsten rods. 
These rods were cut to approximately four centimeter lengths and one end 
of each rod was spot welded to a 0.05 by 0.3 em nickel interface. A 
solid copper lead was spot welded to the other end of the nickel interface. 


This is depicted in Figure 68. 


108 





Tungsten Rod. „Copper Wire 
Nickel 


Interface 


Figure 68. Tungsten electrode connection to the copper lead through 
a nickel interface. 

Two of these electrodes were then cut to lengths differing by one cm 
(2.5 cmand 3.5 cm). It should be noted that tungsten has a tendency to 
split when snipped. The most successful cutting method involves notching 
the tungsten rod with a file and then breaking it over a sharp corner 
mich as glass rod or tubing is cut. The tips were then etched in a 
saturated solution of sodium hydroxide at the same time. The etching 
process is depicted in Figure 69. 


Carbon + AC Power 
Electrode — Supply 


Tungsten 
Probes 





Sodium Hydroxide 


Saturated Solution 





Figure 69. Etching the tungsten electrode tips. 


The etching time was ten minutes and produced identical tips as depicted 


in Figure 70. 


109 





ip diameter .014 CM 





"base diameter 





034 CM 


nA length .1905 cm 


Figure 70. Etched Tungsten Electrode 






In the next step, a five millimeter diameter pyrex tube was cut to 
a length of 20 cm. Two 0.75 mm outside diameter, 0.60 mm inside 
diameter alumina tubes were cut to lengths of 1.5 and 2.5 cm and a piece 
of RG-174/U microcoaxial cable was cut to a length of 30 cm. One end 
of 1E coax was stripped and inserted inside the stainless steel shield 
which is then inserted inside the pyrex tube. At this point the 


electrodes are soldered to the microaxial cable. Figure 71 depicts the 


Rope im this state of completion. 





a m m pa me 






a 


stainless shield 


Figure 71. Nearly Completed Double Probe 


To insure proper insulation, all exposed leads in the two electrode 
circuits are covered with a thin coat of "RIV." When this is dry the 
probe should be checked for circuit continuity. (zero ohms (B) to B') 


and between (A) and (A'). 
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The supporting structure for the elertrode tips is made of liquid 
procelain. The exposed electrode tips are lightly tied together with 
thread and the entire circuit is stretched between two clamps. A thin 
coat of procelain is then applied. After the first coat is co the 
probe may be removed from the clamps. Several additional thin coats 
Should be applied until the desired strength is achieved. Figure 72 
depicts the finished probe. 


LIQUID PORCELAIN COVERING 





| 
a CM» (b) | 
| 

| ON E E | 
| | a] 


Figure 72. Finished Double Probe (a) internal configuration (b) 
external dimensions.. * Insure that the liquid procelain 
is drawn into the aparatus here to provide adequate 
sprenpgch. 

"Because of the larger size of this probe, there was some concem 
that the plasma might be perturbed sufficiently to distort the data. 

To determine what effect this larger probe had on the plasma features a 

comparison was made between the probe response of the finest probe 

available and this slightly larger probe. The comparison is duplicated 
as Figure 73. As can be seen, all features are clearly retained, the 
primary difference being a larger signal from the probe with the larger 


electrodes. Because of the additional size of this probe it was not 


ee) 








used in any experiments where additional diagnostic information was 


taken downstream. 


(a) 
(b) 


—probe response —» 


20 NSEC/500 MV PER DIV. 


Figure 73. Comparison of the density profile sampled at (0, 0, 1.0) by 
(a) the large probe biased at -20 volts and (b) by a very 
fine probe biased at -15 volts. 
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APPENDIX B 
I. PRESSURE PROBE THEORY 


When a force is exerted on a piezoelectric crystal such as quartz, 
it results in the formation of a stress within the crystal. This stress 
distorts the ionic structure of the crystal resulting in an increase in 


the polarization of the sample. This process is depicted in Figure 74. 


STRESS 





STRESS 


(a) —(b) 


Figure 74. (a) Unstressed Piezoelectric Crystal 
(b) Stressed Piezoelectric Crystal 
The stress changes the polarization by AP, the induced 
piezoelectric polarization. 


Recall that the electric displacement vector is given by: 


— z T gp 
(30) D, E P_ 


For a one dimensional strain, the piezoelectric equation is: 


NR = Mer 
(3) D, = Ta, 


When D, = electric displacement vector, T = stress, ds, is the 


Piezoelectric Strain Constant, E = electric field and € is the 


permittivity. 
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For an X-cut quartz crystal [Ref. 18] 


2 725X 10712 coulombs/newton 


(a) à, 


(o) € =4.06x 101 farad/meter 


In the static case the open-circuit voltage is given by [Ref. 4]: 
(32) g=d,/e = 0.055 IF/tw volt-meters/newton 
and the short-circuited charge for a given applied force is: 
(33) Q= dy FP/t 


where F is the force which is considered positive for an extensional 
stress, f = crystal length, w = crystal width and t = the crystal 
thickness. 

The material presented above is satisfactory for low frequency impulses. 
In this investigation, the anlysis is further complicated because the 
plasma produces a stress wave with a short wave length. Graham, et al., 
[Ref. 18] have analyzed the short-circuit piezoelectric current induced 
by the application of a rapidly changing impulsive load. The load 
propagates as a stress wave along the x-axis of a quartz disk. ‘The 


displacement current generated by this stress wave will be: 
(34) i = A aD/at 


Where (A) is the electroded area and D is the displacement vector. 


Substituting equation (31) for (D) we find that: 


f £ f 
(35) far = fio dx + ef Ex dx 
0 0 0 
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Where f is the thickness of the quartz disk. By assuming that a short 
circuit exists between the electrodes and that the crystal has constant 


permitivity the last integral in equation (35) becomes zero [Ref. 18]. 


f 
(36) [oa = 0 


0 
The conductivity of the crystal is essentially zero so that equation 


(35) becomes: 
Po — 
Baal. P(x)dx = D 
fo 
At this point it is assumed that the piezoelectric polarization, 
P(x), is proportional to the x-component of the stress by a coefficient 


(f). This coefficient is independent of time and stress for a given 


Seress range. Iherefore: 
(38) P(X) = fo (x) 


Where c(x) is the x-component of the stress. This assumption is 
verified by Graham [Ref. 19]. With the assumption that the quartz is 


linearly elastic we can write: 
(39) o(x,t) = o(x- Ut) 


Where Us is the wave propagation speed in the crystal. 
The solution for the current from equations (34) and (37) is there- 


fore found to be: 


so 


? o 9X Bb 





/ 
(40) i = adD/at = tau, fao(x)dx = faule, - y ) 


115 





Where 6, is the x-component of stress at the stress input electrode 


and o is the x-component of stress at the rear electrode. 


For early times in the wave transit, i = 0 and equation (40) 
becomes: 
CE SO rll q 


D ——s"0 


L 


Equation (41) a that for times less than the wave transit time, 
the current is proportional to the stress which is proportional to the 
force applied to the crystal which in the present investigation equals 
the momentum transported to the crystal by the plasma. Therefore in 
this investigation the pressure probe response is directly proportional 


to the momentum transport in the plasma. 
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II. PRESSURE PROBE CONSTRUCTION 


Stanford Research Institute provided the pressure probe used in 
this investigation. It is a Sandia Quartz Gauge. The time resolution 
of this type gauge is on the order of a few nanoseconds [Ref. 18]. 


Figure 75 depicts the construction of the gauge. 
X CUT QUARTZ OLD PLATE 












MICRO COAX 


ALUMINUM HOUSING 


Figure 75. Sandia Quartz Gauge. 


This gauge used a guard ring configuration described in detail by 


Graham [Ref. 18]. The configuration is pictured in Figure 76. 
GOLD PLATE 
GUARD RING fi 


GAUGE FACE GAUGE EDGE 





Figure 76. Guard Ring Electrode Configuration 


The inner porticn of the quartz disc is isolated electrically from the 
outer portion by separating the vapor coated electrode into two regions 


as shown above. 
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This configuration helps to eliminate field distortions. At the 
outer edge of the disk, the discontinuity in electric potential and 
dielectric permittivity causes electric field fringing similar to that 
found in parallel plate capacitors. The guard ring eliminates this 
distortion by restricting observation to the central region of the disk 
where there is negligible fringing. 

The guard ring also helps to eliminate distortions due to the 
generation of unloading waves. As the stress pulse moves through the 
disk in the axial direction, boundary conditions generate shear and 
dilatational waves immediately behind the wavefront at the lateral edge 
of the disk. These unloading waves then propagate laterally inward 
from the lateral boundary and the crystal is no longer in a state of 
one dimensional strain. The central region of the disk will be in one 
dimensional strain only for the first wave transit time if the width of 
the outer electrode is such that an unloading wave does not reach the 


central region during the first wave transit time. 
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